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■ Abstract 

(f) ' In this paper we discuss the concept of cosymmetries and co-recursion operators for difference 

t-H . equations and present a co-recursion operator for the Viallet equation. We also discover a new 

type of factorisation for the recursion operators of difference equations. This factorisation 
KH ■ enables us to give an elegant proof that the recursion operator given in arXiv:1004-5346 is 

indeed a recursion operator for the Viallet equation. Moreover, we show that this operator 
c"| . is Nijenhuis and thus generates infinitely many commuting local symmetries. This recursion 

' operator and its factorisation into Hamiltonian and symplectic operators can be applied to 

Ch . Yamilov's discretisation of the Krichever-Novikov equation. 

■ 1 Introduction 

cn : 

The existence of hierarchies of infinitely many symmetries and/or conservation laws is one of the 

■ important characteristics of integrable equations. It has been successfully used in the classification 
\ of given families of partial differential and differential-difference equations [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. 

The hierarchies of infinite symmetries can often be generated by recursion operators map- 
ping symmetries to symmetries [11, 12, 13, 10]. Most of the known examples of recursion opera- 
tors have the property of being Nijenhuis operators. This property was independently studied by 
' Fuchssteiner-Fokas [14, 15, 16] and Magri [17], who named such operators as hereditary symmetries. 

^ It guarantees that one can generate infinitely many local commuting symmetry flows starting from 

some seed symmetries. Furthermore, the conjugated recursion operators can generate infinitely 
many local covariants (the variational derivatives of conserved densities) . For Hamiltonian systems 
recursion operators can be factorised into symplectic and Hamiltonian operators. 
The most well known example is the famous Korteweg-de Vries (KdV) equation 

u t = K = u xxx + 6 u u x , 

which possesses a recursion operator 

sft = D 2 X + Au + 2u x D- 1 , 

where D^ 1 stands for the inverse of D x . The operator ?ft satisfying 

D Ut -K^ = %lD Ut -K , 
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where D ut _K is the Frechet derivative along the KdV, is a Nijenhuis operator and generates the 
local commuting KdV hierarchy 

u tj = W(u x ), j = 0,1,2,- ■■ , 

the first member of which, i.e. a seed symmetry, is ut = u x . 

The operator conjugated to the recursion operator 5ft for the KdV equation is 

3T = D 2 X + Au - 2D~ l u x 

and satisfies the operator equation 

K-K®* = M*D* Ut _ K . 

It recursively produces infinitely many local covariants G n = 5R*(G n _i), n € N, with Go = 1, 
which are also cosymmetries of the KdV equation. Indeed, they satisfy D^ t _ K (G n ) = 0. Thus, for 
the KdV equation the conjugate operator 5ft* is a co-recursion operator mapping cosymmetries to 
cosymmetries. 

The recursion operator 5ft can be represented as the composition of a symplectic operator X and 
a Hamiltonian operator Tii, i.e. 

» = Wi-X, Z = D-\ Hi = Dl + AuD x + 2u x . 

The KdV equation is a multi-Hamiltonian system with infinitely many Hamiltonian operators 
7i n = ^ft n • Ho, n £ N, where Hq = X v = D x . It is not difficult to show that all Hamiltonian 
operators are weakly non-local. The concept of weak non-locality was discussed in detail in [18]. 

In our recent work [19], we adapted the concept of recursion operators for difference equations. 
In particular, we claimed that the operator 9t given in Theorem 5 of [19] (see also (25) below), is 
a recursion operator for the Viallet equation [20]: 

Q := aiit ,o«i,o«o,i«i,i 

+a 2 (ito,o«i,o«o,i + «i,o«o,i«i,i + «o,i«i,i«o,o + «i,i«o,o«i,o) 
+123(^0,0^1,0 + uo,iui,i) + a4(Mi,o«o,i + m ,o«i,i) 

+05(^0,0^0,1 + ui )0 «i,i) + a 6 (u ,o + ^1,0 + «o,i + «i,i) + a 7 = , (1) 

where are free complex parameters. Equation (1) is a quite general and useful example of 
an integrable difference equation. By a point fractional-linear transformation this equation with 
a generic choice of parameters can be reduced to Adler's equation, also referred to as the Q4 
equation in the ABS classification [21, 22]. All of the ABS equations can be obtained from the 
Viallet equation by a simple specialisation of the parameters. Thus, the results obtained for the 
Viallet equation can be easily applied to all ABS equations. 

In this paper, we introduce and discuss the concepts of cosymmetries and co-recursion operators 
for difference equations. We give a complete proof of Theorem 5 formulated in [19], namely that 
the operator 9^ presented in the theorem is indeed a recursion operator of the Viallet equation (1) 
and it is a composition of Hamiltonian and symplectic operators. Moreover, this operator possesses 
the Nijenhuis property and thus produces infinitely many commuting local symmetries starting 
from two seed symmetries. Equation (1) possesses an infinite hierarchy of cosymmetries which can 
be generated by a co-recursion operator. Here we would like to stress that the cosymmetries of 
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equation (1) do not coincide with its covariants and that the co-recursion operator is not conjugated 
to the recursion one. As a by-product, we have found a recursion operator and its factorisation into 
Hamiltonian and symplectic operators for the differential-difference integrable equation known as 
Yamilov's discretisation of the Krichever-Novikov equation [10]. 

In order to be self-contained we introduce our notations, sketch the framework and give ba- 
sic definitions in Section 2. More detailed description of the framework and motivations for the 
definitions the reader can find in our previous paper [19]. In Section 2.1 we define the field of 
fractions Tq corresponding to a difference equation Q = 0, the elimination map £ and the set of 
dynamical variables. In the next subsection we give the definitions of symmetries and cosymmetries 
illustrating by useful examples for equation (1). 

In Section 3 we discuss difference and pseudo-difference operators. In particular we define 
weakly non-local pseudo-difference operators. They are difference analogues of weakly non-local 
pseudo-differential operators in the theory of Poisson brackets (see A.Ya. Maltsev and S.P. Novikov 
[18]). We give the definitions of recursion and co-recursion operators for difference equations and 
show that operator £K (25) is a recursion operator for the Viallet equation (1). On the course of 
the proof we discover a new factorisation property of the recursion operator (25). Namely, there is 
a constant jjl (33) such that the operator 9\ — fi can be factorised over the field Tq into a difference 
and a weakly non-local pseudo-difference operator. Interestingly enough, the coefficients of £K are 
elements of the subfield T s C Tq (of rational functions of variables {ti n ,o I n £ see Section 
2.1) while the coefficients of the factors are in Tq and cannot be restricted to T s - A co-recursion 
operator for equation (1) is given explicitly in Theorem 2. 

In Section 4 we discuss Hamiltonian and symplectic operators related to difference equations and 
show that the recursion operator (25) is hereditary (Nijenhuis). The main result of this section is 
Theorem 3, which states that symmetries of equation (1) generated by the recursion operator £K (25) 
are all local and commuting, despite of d\ being a weakly non-local pseudo-difference operator. The 
definitions of Hamiltonian and symplectic operators, the hereditary property of a recursion operator 
as well as the proofs of assertions are not different from the continuous case if we reformulate the 
problem in terms of a variational complex. It enables us easily to adapt the theory developed for 
the continuous case [23, 24, 7, 25] to the difference one. 

2 Framework and basic definitions 

To make the paper self-contained, we give the basic definitions of the elimination map, dynamical 
variables, symmetries and cosymmetries of difference equations. Motivations for these definitions 
have been discussed in detail in our paper [19]. 

2.1 Difference equations and dynamical variables 

Difference equations on Z 2 can be regarded as a difference analogous of partial differential equa- 
tions. Let us denote by u = u(n, m) a complex-valued function u : 1? i-> C, where n and m are 
"independent variables" and u will play the role of a "dependent" variable in a difference equation. 
Instead of partial derivatives we have two commuting shift maps S and T defined as 

S : u i-> u±fi = u{n + 1, m), 
T : mi-> uo,i = u(n, m + 1), 
S p T q : u i-> u ps = u(n + p,m + q). 
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For uniformity of the notation, it is convenient to denote the "unshifted" function u as u 0j o- In the 
theory of difference equations we shall treat u Ptq as variables. The set of all shifts of the variable u 
will be denoted by U = {u p>q \ (p, q) G Z 2 }. For a function / = f(u pijqi , . . . , u Pkm ) of variables u p>q 
the action of the operators S, T is defined as 

S % T ] (f) = fij = f (u Pl +i tqi +j , . . . ,U Pk +i,q k +j). 

In this paper, we consider a quadrilateral difference equation of the form 

Q{UQ,0, Ulfl, «0,1) = i ( 2 ) 

where Q(uofi, ui ; o, uo,i, ui,i) is an irreducible polynomial of the "dependent variable" u = uo,o and 
its shifts. Irreducibility means that Q cannot be factorised and presented as a product of two 
polynomials. It is assumed that equation (2) is valid at every point (n, m) G Z 2 and thus (2) 
represents the infinite set of equations 

Qp,q = Q(Up,q,U p +l,q,U p ,q+l,Up+l,q+l) = 0, (p, q) G Z 2 . (3) 

Moreover, we assume that Q is an irreducible amne-linear polynomial which depends non-trivially 
on all variables, i.e. 

d Ui jQ^0, dl..Q = 0, i,je{0,l}, QeC[Mo,o,Mi,o,Mo,i,Mi,i]- (4) 

Here and later in the text d Ui j f refers to the partial derivative of / with respect to Uij which we 
also often denote as f Uij . 

Let C[U] be the ring of polynomials of the variables U. Maps S and T are automorphisms 
of C[U] and thus C[U] is a difference ring. We denote Jq = ({Q p , q }) the ideal generated by the 
difference equation and all its shifts (3). It is a prime difference ideal and thus C[U]/Jq is a 
difference quotient ring without zero divisors. The corresponding field of fractions we denote as 

T Q = {[a]/[b]\a,b€C[U], b £ J Q } , 

where [a] denotes the class of equivalent polynomials (two polynomials /, g G C[U] are equivalent 
if / — g G Jq). For a, b, c, d G C[f7], 6, d Jq, the two rational functions a/b and c/d represent the 
same element of J-q if ad — be G Jq. The field J"q is the main object in our theory. 
The fields of rational functions of variables 

^ s = {u n ,o I n G Z}, C/t = {n ,„ | n G Z}, U = U s UU t . 

are denoted respectively as 

F B = C(U S ), F t = C{U t ), J~o = C(U ) . 

Since the ideal Jq is generated by affine- linear polynomials Q p , q , we can define the elimina- 
tion map 6 which simplifies computations modulo the ideal drastically. One can uniquely resolve 
equation Q = with respect to each variable 

«o,o = F{u-ifi, uo,i, «i,o = G(no,o, «o,i> u i,i)> /g\ 

u ,i = ^(«o,o,«i,o,«i,i), ui,i = M(uo fi ,ui )0 ,u 0) i) 

and since Q is an affine linear polynomial, functions F, G, H and M are rational functions of their 
arguments. Equations (5) enable us recursively and uniquely to express any variable u ps in terms 
of the variables Uq = U s U Ut . 
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Definition 1. For elements of U the elimination map £ : U i-> C(Uq) is defined recursively: 

Vp€Z, £ (u ,p) = u , p , £(« P ,o) = « P ,o , 

if p > 0,g > 0, £ (u Pj9 ) = M(£ (u p -i :q -i),£(u P:q -i),£(u p -i, q )) , 

ifp<0,g>0, £ («p j9 ) = H(£ (up j g_i),£(u p+ i ) g_i),f (up+i,,)) , (6) 

if p > 0, q < 0, f (« P , g ) = G(£(up-i jq ),£(up-i >q +i),£(up !q +i)) , 

if p < 0, q < 0, £(u P ,q) = F(£(u p+ i, q ), £(u p , q+ i),£(u p +i, q+ i)) . 

For polynomials f(u Pim , . . . ,u PktQk ) € C[U] the elimination map £ : C[U] i-> C(J7o) ^ s defined as 

£ : f(u pi 

,qii ■ ■ ■ ' U Pk,qk \ u Pk,qk )) G C(C/ ). 

For rational functions a/b, a,b G C[f7], b Jq the elimination map £ is defined as 

£ : a/b^ £(a)/£(b). 

It follows from the above definition that £ : Jq *— > 0. The elimination map is the difference field 
isomorphism £ : Tq i— > J-"o- Any element of J-q can be uniquely represented by a rational function 
of variables Uq. Variables Uq we shall call the dynamical variables. 

2.2 Symmetries and cosymmetries of difference equations 

Symmetries and conservation laws of difference equations have been discussed in detail in our paper 
[19]. In this section we recall the definition of a continuous symmetry of a difference equation 
and give the definition of a cosymmetry. For difference equations cosymmetries do not coincide 
with covariances (the variational derivatives of conserved densities) as in the case of evolutionary 
equations. Cosymmetries are new objects in the theory of difference equations, and they will play 
important role in our construction of recursion and co-recursion operators. 

Definition 2. Let Q = be a difference equation. Then K € Tq is called a generator of an 
infinitesimal symmetry ( or simply, a symmetry ) of the difference equation if 

Dq(K) = 0. (7) 

Here Dq is the Frechet derivative of Q defined as 

Dq=Y, Qu^S'V , Q UiJ = d Ui .Q. (8) 

For a quadrilateral equation the sum in (8) has only four terms. In (7) Dq(K) is equal to zero 
as an element of Tq. The way to check this is to apply the elimination map and therefore K is a 
symmetry if £(Dq(K)) = 0. If the difference equation Q = admits symmetries, then they form 
a Lie algebra denoted as 21q [19]. 

Similarly, we define cosymmetries: 

Definition 3. Let Q = be a difference equation. Then uj G Tq is called a cosymmetry (or a 
characteristic of a conservation law ) of the difference equation if 

D* Q (u) = 0, (9) 

where Dq is a formally adjoined operator of the Frechet derivative which is defined as 

D* Q = ^ 5 "T '(}„ ,. (10) 
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Cosymmetries of equation Q = form a linear space over C which we denote as Oq. 
We shall see that equation (1) has infinite hierarchies of symmetries and cosymmetries which 
can be naturally expressed in terms of the discriminants of Q 

h(uo j0 ,uifl) := Qd UQ1 d Ull Q - d UQ1 Qd Ull Q; (11) 

9(ui,o,u ,i) ■= Q d U0 d ul l Q - d U0 Q d Ul l Q ; (12) 

and the function 

w = . (13) 

For any / G Tq we shall use the abbreviated notation f k = S k (f). In particular 

h k = S k h(u 0fi ,u lfi ), g k = S k g(u lfi , u ,i), w k = S k w, 

In this notation, the first two symmetries of equation (1) (see [26, 27]) can be written in the 
following form 

K {1) =hw- -d ui0 h = h-iw + ^-1,0^-1 ; (14) 
K {2) = h h^w 2 ( Wl + w-i) . (15) 
In [28, 27] it was shown that the above symmetries commute 

[K^,K^] = D Km (KW) - D KW {K<®) = 0. (16) 

It is much less known about cosymmetries. For partial differential equations cosymmetries were 
introduced and studied in [7]. The simplest cosymmetry for equation (1), which we have found by 
direct solution of equation (9), can be written in the form 

w (i) = wi _ w + d ui0 g + d UQQ h 

It is easy to check that a/ 1 ) is a cosymmetry. Indeed, using the elimination map one can show that 

K-(i) K W 

s-^q^J^ + uWq^ = -t-'s-HQu^) - r-V 1 ^) = -±- - J^- (18) 

This immediately leads to DqOJ^ = and proves that a/ 1 ) is a cosymmetry of (1) according to 
Definition 3. The next cosymmetry oj^ (cf. (39)) looks considerably more complicated and it will 
be derived in the next Section. 

Symmetries are elements of F s , they do not depend on variables u k)V with p 7^ and 

thus it is relatively easy to find them starting from the Definition 2. Cosymmetries a/ 1 -*, oj^ are 
elements of J-q which cannot be reduced to elements of J- s or Tt by the elimination map. 

Equation (1) has an infinite hierarchy of commuting symmetries as well as an infinite hierarchy 
of cosymmetries. In the next sections we shall show that the hierarchy of symmetries can be 
generated by successive applications of a recursion operator to the seed symmetries K^ 2 \ while 
the hierarchy of cosymmetries can be obtained by applying successively a co-recursion operator to 
the seed cosymmetries 
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3 Recursion and co— recursion operators of difference equations 



In this section we prove that equation (1) possesses weakly nonlocal pseudo-difference recursion 
and co-recursion operators. 

We remind that difference operators are finite sums of the form 

A = a n S n + a^S"" 1 + ■■■ + a m S m , a k G Tq . (19) 

They form a non-commutative ring J^f^S] of polynomials in S with coefficients from Tq. Difference 
operators we also call local operators since the action A : a t- > A(a) G Tq of a difference operator 
A is naturally defined for any element a G Tq. 

For our theory the operator of a finite difference A = S — 1 is of particular importance. This 
is a difference analogue of the derivation D x in a differential field in the continuous case. The 
corresponding inverse pseudo-difference operator A^ 1 is the analogue of D~ x or integration. The 
action of the operator A -1 is not defined for all elements of Tq, but only on the elements of 
image space ImA = ATq. In general the kernel space Ker A C Tq could be rather nontrivial (see 
discussion in [19]), but in the case of equation (1), where Q is an irreducible polynomial satisfying 
(4), one can show that Ker A = C. Thus for a = b\ — b, be Tq we have A -1 (a) = b + a, a£C. In 
what follows we shall ignore "the constant of integration" a, since its effect on results is inessential 
similar to the continuous case [2]. 

We define weakly nonlocal pseudo-difference operators as finite sums of the form 

5 = 5o+a 1 A- 1 o& 1 + a 2 A- 1 o&2 + --- + a m A- 1 o6 m , a k , b k G Tq, B G Tq[S\. (20) 

It is a difference analogue of weakly nonlocal pseudo-differential operators which play important 
role in the theory of multi-Hamiltonian partial differential equations [18]. A composition of two 
weakly nonlocal pseudo-differential operators is a pseudo-difference operator, but not necessarily 
weakly nonlocal. 

Formally conjugated difference (19) and weakly-non local pseudo-difference (20) operators are 
defined as 

A* = S~ n o a n + S 1 - 11 o a n _! + • • • + S~ m o a m , 
B* = Bq — hA^S ooi- b 2 A- 1 S oa 2 bmA^S o a m . 

The action of a pseudo-difference operator is not defined for all elements of Tq. The subset of 
Tq for which the action of B is defined Dom(B) = {a G Tq \ B(a) G Tq} is called the domain of 
B. For example Dom(A _1 ) = ImA and for A G Tq[S] we have Dom(A) = Tq. 

By a recursion operator of a difference equation (1) we shall understand a pseudo-difference 
operator 9^ such that 

^:Dom(^H)n2l Q ^2lQ, 

where 21q is the linear space of symmetries of this difference equation. In other words, if the action 
of 91 on a symmetry K G Tq is defined, i.e. 9K(K) G Tq, then 9\(K) is a symmetry of the same 
difference equation. The operator of multiplication by a constant is a trivial recursion operator. 
It easy to see that a pseudo-difference operator %\ is a recursion operator for a difference equation 
Q = if there exists a pseudo-difference operator such that 9\ and *P satisfy the following 
operator equation 

D Q o<R = <#oD Q . (21) 
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Indeed, if K is a symmetry of this difference equation and K = 9t(K) G Tq, then it follows from 
(21) that DqK = 0, and thus K is also a symmetry. 

The equality in (21) has to be understood in the sense of the field J-q, or in the usual sense 
after application of the elimination map to the coefficients of the pseudo-difference operators. In 
[19] it was shown that the operator equation (21) is equivalent to two equations 

o (Q U1Q S + Q UOi0 ) = (Q U10 S + Q U0fi ) o <K , (22) 

and 

V o (Qu hl S + Q U0A ) = {Q Uia S + Q U0A ) o T(9t) . (23) 

We say a pseudo-difference operator 2U is a co-recursion operator of a difference equation 
Q = if it maps cosymmetries to cosymmetries 2U : Dq nDom2U h-> Oq. If the difference equation 
possesses a recursion operator satisfying equation (21) then 2U = Indeed, conjugating (21) 
we get 

Dq o *p* = 5H* o D* q (24) 

and thus if a; is a cosymmetry and Cj = 2B(cj) € Dom2B, then Dq(6j) = 0. 

In the rest of this section, we show that with the operator 9\ presented in paper [19] there exists 
(and will be given explicitly) a pseudo-difference operator satisfying equation (21). It implies 
that 9\ is a recursion operator of the difference equation (1) and 2U = *p* is a co-recursion operator. 

We rewrite operator 9\ in [19] in the form 



h h-x w 2 wjS 2 + h h-i w 2 w 2 _ 1 S- 2 + 2K^K^ (^S + S' 1 ^ 
-w 2 (h-! h lW 2 + h^hwl,) + (K^K^l + K^K^ 

^^A-'otiA L^A-'otiA (25 ) 

\ h-i hi \ h-i h J 

where h, w, and are given by (11), (13), (14) and (15) respectively. It is a weakly 

non-local pseudo-difference operator and can be written as the product of the following operators 



h-\hh\ i h^\hh\ 



(u 1>0 - m-i,o) 2 (^2,o - UQ fl ) 2 ' (ui fi - u_i i0 ) 2 (n 2 ,o - «o,o) 2 

+ 2 A _1 S o id 2 ) + 2 A' 1 o K<U , (26) 
I= I 5 -5-'i, (27) 

that is, 91 = Hoi. In section 4, we are going to show that operator H is Hamiltonian and operator 
1 is symplectic. 

Theorem 1. Operator 9\ given by (25) is a recursion operator of equation (1). 

Before we prove this theorem, we first give a few lemmas. 
Lemma 1. If a pseudo- difference operator 

A = a^S 2 + aP>S + + a^S' 1 + a^S~ 2 + b^A'%^ + b {2) A~ l b^ 
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can be factorised as the product of the following two operators 

C = cWs + c^ + c^S^ + c^S^ + b^A^d^ + b^A-'d^, 
E = eWs + eW 

that is, A = C o E, then the coefficients of the operator C satisfy the following conditions 

<F>=<F>£\ a<- 2 > = cC- 2 >ei3, (5 

a (l) = c (0) e (l) + c (l) e (0) j a (-l) = c (-2) e (l) + c (-l) e (0) { r 

6( 2 ) = (4 2 ) e «)_i + 4 2 ) e w ) (j 

a (0) = c (-i) e (i) + 6 ( 1 )(4i) e (i) ) _ 1 + ^(d^eW).! + c (0) e (0). (J 

Proof. Expanding the composition of operators C and E and using the identity A _1 <S = 1+A 
we get: 

CoE = cWeJV + fc(°)e« + c^e^) 5 + ( C (- 2 ) e « + c*" 1 ^ + c^e^" 2 



+ A- 1 ((^eW).! + 4 1 > e <°>) + ftp* A" 1 ((4 2 )e«)-i + 4 2 >e(°>) • 

Comparing the coefficients of the operator A with the above operator we obtain (28)— (32). 
Lemma 2. Let Q and 6e given by (1) and (25) respectively. Then for the constant 

H = (04 — a3)(a3aia7 — 03a 2 — 2a7a 2 + 4asa2a6 — 2aiag + 040107 — 04a 2 ) , (c 
the operator £R — /j can 6e factorised over Fq as 

*R- M = X- (Q Uli0 «S + Quo, ), (J 

w/iere 

M = c (1) 5 + c(°) + c^S- 1 + c(- 2 )«S- 2 + 2 K« A" V 2 > - 2 IfP) A" V 1 ) (c 
iwit/i coefficients 

^ _ hh-iw 2 wf _ (_ 2 ) _ /i /i-i u> 2 w 2 ! _ ^ r 
<^(Qui,o) ^~ 2 (Quo,o) 

(0) _L_ / 2ifWjf( 2 ) _ hh-i w 2 w 2 S(Q U0fi ) \ 

Qu lfi \ h S(Q U1 J J ' ^ 

1 (2KMkW hh^ l w 2 w 2 _ 1 S- 2 {Q Ul0 )\ 



if( 2 ) kJ 2) i<f } 5(Q UOj0 ) K[ 2) uW k[ 2) kW 



(J 



^ {2) = - , - - + — ■ — — + 



hQ ul}0 hQ U00 hiQ ul S(Q U10 ) K W hK[ 1] Q Ul 
H h x w 2 w{ ( hS^jQ^J h_i \ 



2K[ 1) Q u1>0 2K[ 1) \Qu , S 1 (Qu () ,o) Qui,o 



hwZw 2 ( hiSiQnn n) S 2 (Q,,, n „1 



and oj^ given by (17). Moreover, u/ 2 ) is a cosymmetry of equation (1). 

Proof. Using Lemma 1, we can easily determine cW,c(°),c( _1 ) and c*~ 2 ) given by (36-38). It 
follows from (18) that a/ 1 ) (17) satisfies relation (30). Using condition (32), we get 

= f- w j ( h h2W l + h _, hl w 2) + ^ ^(D K (2) + K f) K ^ _ ^ 



1 «S(Qu .o) / 2if{ 1) *f ) _ M«?^\ 
2if 1 (1) Q Uli0 5(Q„ 1 ,o) V >H <S 2 (Q U1 , ) J' 

which can be rewritten as in (39). Substituting the latter into condition (31) 

K (2) K (2) 

we obtain the value of fi. Using the elimination map and Definition 3 one can directly check that 
<jj( 2 ) given by (39) is a cosymmetry of equation (1). ■ 

Factorisations (34) is quite remarkable. Indeed, the coefficients of the weakly nonlocal pseudo- 
difference operator *H — fi are in the field T s , but the coefficients of both of its factors depend on 
the variable 1x0,1 £ F s . It is indeed a factorisation over Tq. 

Similarly, we obtain the following factorisation for operator T(9 c t) — fi: 

Lemma 3. For operator D\ given by (25)and the constant fi given by (33), the following factorisa- 
tion holds for all solutions of equation (1) 

T{V\)-ix = M ■ (Q UM <S + Quo tl ) , (40) 

where 

M = c&S + c<°> + c^S- 1 + c^S~ 2 - 2 kWA-'JV + 2 ^A" 1 ^ 1 ). (41) 
Here uj^ and u/ 2 ) are given by (17) and (39) respectively, and 

= T(hh^w 2 wj) ^ g( _ 2) = T(hh^ 1 w 2 w 2 _ 1 ) ^ ^ 

<0) _ J_ / 2T(KWK( 2 )) _ T(hh^w 2 w 2 )S(Q U0:1 ) \ 

" Q« M ^ T(h) S(Q Uhl ) J' { ] 

.(-D = 1 ( 2T(K^KW) _ T(hh^w 2 w \)S- 2 {Q Ul:i ) \ 

~ s-i(Q U0A ) { r(/i_i) s- 2 (Q U(hl ) )■ 1 ) 

Proof of Theorem 1. Using the elimination map we can directly show the operator identity 

(Qui,o«S + Qu 0}0 ) M = (Qu^S + Q U0A ) M. 
We denote the above operator as From Lemmas 2 and 3, it follows that 

*P ° (Q Ulfi S + Q U0fi ) = {Qu lfi S + Q U0fi ) o (JR - n). 
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and 

*P o (Q„ M S + Q U0A ) = (Q U1>1 S + Q UOjl ) o (T(5H) - /*)■ 

Thus we proved formulas (22) and (23) for the pseudo-difference operators 9^ — /i and *p. Therefore, 
operator !EH — /U is a recursion operator of equation (1). So is operator JH since is a constant. ■ 
From the proof of Theorem 1, it follows that 

D Q o(K-n)=yoD Q (45) 
for all solutions of equation (1). This leads to the following result: 
Theorem 2. Operator 2B = is a co-recursion operator of equation (1), where 

¥ = (Q Uh0 S + Q UOfi ) M 
and operator A4 is defined in Lemma 2. 



4 Locality of symmetry hierarchies 

In this section, we prove that recursion operator (25) is a Nijenhuis operator. Despite of being 
weakly non-local, it generates an infinite hierarchy of local commuting symmetry flows. The main 
result of this section is a re-statement of the Theorem proven in the continuous case [25] with the 
adjustments to the difference variational complex. 



4.1 Difference variational complex and Lie derivatives 

Notice that all coefficients of operator are elements in F s . In this section we sketch the difference 
variational complex over the ring T s [29, 30] in the same spirit as the variational complex [23, 24, 
31, 7]. We also adapt the formula of the Lie derivatives along evolutionary vector fields. 

Recall that T s is the field of rational functions of variables U s = {u n fi \ n G Z}. Without causing 
any confusion, we denote u n $ by it n for simplicity in what follows. 

The field F s is also a linear space over C. Let us consider the extended linear space C s = 
-^s © Span c {log a | a € J-'g}, so that elements of C s are linear combinations with complex coefficients 
of elements in J- s and logarithms of elements in . For any element 5 £ £ s ,we define an equivalence 
class (or a functional) fg by saying that two elements g, h € C s are equivalent if g — h G A(£ s ). 
The space of functional will be denoted by F' s , it is a linear space over C and it does not inherit 
a ring or field structure of J- s . 

The evolutionary vector fields over the ring T s 



du r , 



form a Lie algebra denoted by rj. The action of any element P G I) on fg G J 7 ^ can be defined as 

p°Jg = J xp(g) = I E( 5np )|^ = / D 9\n m 



This is a representation of the Lie algebra f). We build up a Lie algebra complex associated with 
it. This complex is called the difference variational complex. Here we give the first few steps. 
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We denote the space of functional n- forms by Sl n starting with f2° = T' s . We now consider 
the space Q 1 . For any vertical 1-form on the ring P s , i.e., uj = J2 n h^du n , there is a natural 
non-degenerate pairing with an element P G f): 

<u, P>= j ^h^S n P = j I^S^h^A P . (47) 

Thus any element of fi 1 is completely defined by £ duo = ^2 n S~ n h^duo. We simply say £ G fi 1 . 

The pairing (47) allows us to give the definition of (formal) adjoint operators to linear (pseudo-) 
difference operators [29, 31]. 

Definition 4. Given a linear operator A : \) ^ ^l 1 , we call the operator A* : f) — > the adjoint 
operator of A if < APi, Pi >=< A*P 2 , Pi >, where Pi G f) for i = 1, 2. 

Similarly, we can define the adjoint operator for an operator mapping from Q 1 to f), from f) to 
f) or from to fi 1 . 

The difference variational derivative (Euler operator) of each functional jg G P' s denoted by 
6 U0 (jg) G SI 1 is defined so that 

< S U0 ( j g), P>=(dj 9)(P) =P° Jg = J ( D g [P] = j D*(1)P 

•&)'-<E*-&.'>. 

/ nSZ 

where d : £l n — > Q n+1 is a coboundary operator. Due to the non-degeneracy of the pairing (47), we 
have 

For any £ G fi 1 , it can be shown that d£ = — [29, 24]. Thus d£ = is equivalent to 
= D|. We say that the 1-form £ is closed if d£ = 0. 

We can define the Lie derivative of a given object in the complex along vector field K G f). It 
can be expressed explicitly in terms of Frechet derivatives as follows [24]: 

L K Jg = jD g [K] for jg G P' s ; 

L K h = [K,h] for h G f); 

1^ = 1^+1^(6 for (G!) 1 ; 

L^H = Ah[K] - D K m + <RD K for JH : F) -> fj; 1 J 

L^-H = - D^ft - for ft : Q l -»• f); 

L^Z = Djjft] + D^-Z + XDk for X : F) -> fi 1 . 

For a given object a and K G f), if Lft-cx = we say a is conserved (or invariant) along the vector 
field K and the vector field K is a symmetry of the object a. This complex plays an important 
role in the study of differential-difference equations or symmetry flows of difference equations since 
we can associate equations (flows) with vector fields. Hence we can define interesting objects such 
as conserved densities, symmetries and recursion operators as invariants along vector fields. 
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4.2 Symplectic, Hamiltonian and Nijenhuis operators 

In the context of difference variational complex, we now define symplectic and Hamiltonian opera- 
tors. Such definitions are the same as in the continuous case, where we define them in the context 
of the complex of variational calculus. All the results to determine whether a given (pseudo-) dif- 
ferential operator is Hamiltonian or symplectic are still valid. We shall use them later on without 
redeveloping them. 

Definition 5. A linear operator A: f) ^ Q 1 (or Q, 1 — > t)) is anti- symmetric if A = —A*. 

Given an anti-symmetric operator I : t) — > fi 1 , there is an anti-symmetric 2- form associated 
with it. Namely, 

u(P,G)=<X(P),G>=-<X{G),P>=-u{G,P), P,Get). (50) 
Here the functional 2-form oj has the canonical form [31] 



X - J diiQ A IdiiQ. (51) 



U 2 

Definition 6. An anti- symmetric operator Z : f) — >■ Q 1 is called symplectic if its associated 2-form 
(51) is closed, i.e. du = 0. 

Proposition 1. For any function f G P s depending only on uq and u\, the operator fS — S~ l f is 
symplectic. 

Proof. This operator is obviously anti-symmetric and its associated canonical 2-form is 



= \ J (f du o A dn i ~~ duQ A S 1 (fdu )) = J /du A dui . 



Now we have 



// q ^ d f "\ 

I — — duo A duo A du± + — — dni A duo A du\ I = 0. 
\du dux J 

Therefore a; is a closed 2-form and thus the given operator is symplectic. ■ 
It immediately follows that operator I defined by (27) is symplectic since the difference poly- 
nomial h given by (11) depending only on uq and u±. 

For an anti-symmetric operator % : VL 1 — > h, we define a bracket of two functionals as 

{//, fy} =< S U0 (Jf), nS U0 (Jg)>. (52) 

Definition 7. The operator % is Hamiltonian if the bracket defined by (52) is Poisson, that is, 
anti- symmetric and satisfying the Jacobi identity 

{{MM, M+{{fyjh}, m+{{jhjf}, j 9 }=o. 

The Jacobi identity is abstractly the same as in the continuous case when % is a differential 
operator. Thus the results to determine whether a given difference operator is Hamiltonian are still 
valid. In [31] (see p. 443), it is formulated as the vanishing of the functional tri- vector: 

X H{e) (@ n ) = 0, where O n = ± j 6 A H{0) . (53) 

It is clear that any anti-symmetric constant operator is Hamiltonian. For example, operators 
S — S^ 1 , (S + 1)(<S — 1) _1 and (S — 1)(S + 1) _1 are all Hamiltonian operators. 
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Proposition 2. Anti- symmetric operator H defined by (26) is a Hamiltonian operator. 

Proof. According to (53), we check X^(0)(0%) = 0, where the associated bi-vector of 7i is 

&n = \ J 0AH{9) 

= j 9 A (h-! h h^wlOx + 2 A^K^e} . 

Instead of writing out the full calculation, we only demonstrate the method by picking out terms 
of 9 A A^K^O A A- l K^9. The relevant terms in X H{e) (@ n ) are 

f (0 A X n(e) (K^) A A^K^O + K^6 A A-\X n(e) (KW) A 9) 
= J 9 A X n(e) (K^) A A~ 1 K^9 + JiS' 1 - 1)^X^9 A X n(e) (K^) A 9 
= J6A X H(e) (K^) A A- l K^9 - J A- l K^9 A X H{e) (K^) A 9 (54) 
= j 9 A X n(e) (K^) A A- l K^9 + / 9 A X n(e) (K^) A A~ l K^9 
= 1° AEi^(^(^)9 Ui K( 2 ) AA~ 1 K^e + j9 AE i 5 i («(e))a uj K( 1 ) A A- l K^9 . 

Notice that the terms with either A" 1 ^ or A" 1 ^ 2 ^ in S i (H{9)) for i G Z equal to 

&(H{0)) : 2K^ ) A- 1 K^9 + 2k\ 2) A- 1 K^9. 

Substituting it into (54), we obtain that the coefficient of 2 9 A A^K^ 2 ^ A A^K^O is: 

d k(2) [kW] - d kW [kW] = [K^\ k<V\ = . 

By working out for other terms, we can show that the tri- vector X^(#)(@%) vanishes, which implies 
that rl is a Hamiltonian operator. ■ 
The Jacobi identity is a quadratic relation for the operator %. In general, the linear combination 
of two Hamiltonian operators is no longer Hamiltonian. If it is, we say that these two Hamiltonian 
operators form a Hamiltonian pair. Hamiltonian pairs play an important role in the theory of 
integrability. They naturally generate Nijenhuis operators. 

Definition 8. A linear operator d\ : f) — > f) is called a Nijenhuis operator if it satisfies 

[<rp, %\g] - m[mp, g] - <k[p, *hg] + m 2 [p, g\ = o, p, g e r>. (55) 

Using formula (49), this identity is equivalent to 

L^gK = KL G K, Gel). (56) 

An equivalent formulation is: D^fKP\{G) — ^J{D<^[P}(G) is symmetric with respect to P and G 
[16]. It can be used to check directly whether a given operator is Nijenhuis or not [25]. 

The properties of Nijenhuis operators [24] provide us with the explanation of how the infinitely 
many commuting symmetries and conservation laws of integrable equations arise. In applications, 
there are nonlocal terms in Nijenhuis operators. For pseudo-differential operators, a lot of work has 
been done to find sufficient conditions for Nijenhuis operators to produce local objects [32, 33, 25]. 

Proposition 3. Operator 9\ defined by (25) is a Nijenhuis operator. 



14 



Proof. We need to check that the expression H := D<r[9\P](G) — 9KD<r[P](G) is symmetric 
with respect to P and G. The calculation is straightforward but rather long and not suitable for 
a presentation in a journal article. Here we show only one step by picking out the terms in H 
involving either G4 or P4. We use the notation Pj = S l P and Gj = S^G. Since 9\ is a second order 
operator, we compute the terms containing either G2 or P2 (P3 for nonlocal terms) in expression 
Ah[P](G). These terms are 

-2hh-! w 3 wl (Pi - P_i)G 2 - 2hh-! w 2 wf(P 2 - P)G 2 + D h [P]h-! w 2 w\G 2 
+hD h _ 1 [P] w 2 wlG 2 - 2K {1) h- 1 w 2 wlP 2 G 1 - 2K {1) hw 2 wlP 2 G- 1 + 2w 2 h_ 1 h lW \P 2 G 

-h^w 2 w 2 (<9 U2 /ii) P 2 G - 2K a) 1 hw 2 w 2 1 P 2 G + 2 A" 1 (hiw\wlP 3 G) 

-2w 2 hh^ 1 w 2 1 P 2 A' 1 (K^G/h-! - K^G/h] . (57) 



Now the term hh-\ w 2 w 2 S 2 in 9\ acting on the above expression leads to the terms with either G4 
or P 4 in 9\D<r[P](G), denoted by H 1 ' 4 . We obtain 

H 1 ' 4 = -2h_ x hh x h2 w 2 w\w\w\{P 3 - P 1 )G 4 - 2h- 1 hh 1 h 2 w 2 w\wlwl{P A - P 2 )G 4 

+h^ 1 hh 1 w 2 w\w\w\D h2 [P] G 4 + h-ihh 2 w 2 w\wlwlD hl [P] G 4 - 2h^ 1 hh 1 w 2 w 2 l K < 2 l) wlwlP i G 3 

O O ( 1 1 O O O O O Q O O O / *-\ 

—2hh-\h,2W w 1 K 2 w 2 w 3 PiGi + 2hh-\h\h 3 w w 1 w 2 w 3 PiG2 — /i-i /i/iiw u^u^w^ (3 M4 /i3) P4G2 

-2/i/i_iu; 2 ^i ; r{ 1) /i2u;lw|P4G'2 + 2/i /i_i/i 2 w 2 u> 2 u/^fPtGi 

-2/i 2 u> 2 w 2 u^ 2 P 4 (/ii/i_iK (1) Gi - h^hK^d + h^K^G - h^K^G^j 



0000/ x _ \ 00(1)00 

-2h-\hh\ w w 1 w 2 w 3 I I12W2 — -d u:i h2 I P3G4 — 2h^\hh\ w w-^K^ w 2 w 3 P l iG 3 



2 

+2hh-\hi w 2 w\wlw 2 ^h 2 w 3 + ^5 U2 /i 2 ^ P2G4 + 2hh_ 1 h 1 w 2 w\w 2 w 2 ^h 3 w 3 - ^<9 U4 /i 3 ^ PtG 2 

0000/ 1 _ \ oofl^oo 

+2/i_i/i/i2 w w 1 W2W 3 I /iiW2 + —o Ul hi J P1G4 + 2ft./i_i/i2 iu 

Q Q Q Q OOOOl 1 \ 

+h-\hh2W w 1 w 2 w 3 (d U2 hi) P2G4 — 2hh^ih,2W w l w 2 w 3 Ihiwi — -o U2 h\ I P4G2 

-2/i 2 w 2 w 2 ^ 2 P 4 (/ii/i-iK (1) Gi + hhxK^G - h-^K^G^j . 

We then collect the terms with either G4 and P4 in D<r[9\P](G), denoted by H 2 ' 4 . These terms 
can be obtained from the terms containing P2 in (57). Simply replacing P2 by <S 2 (91P), we get 

-2/i A? (9*P) 2 G 2 - 2K ( - 1 h_ 1 w 2 w 2 1 (<RP) 2 G 1 - 2K^hw 2 w 2 (<KP) 2 G- 1 
+2w 2 h^ h lW f(m.P) 2 G - w 2 h^ (a u2 /ii) (<RP) 2 w 2 G - 2K ( L ) 1 hw 2 w 2 (mP) 2 G 
+2 A' 1 (h lW 2 w 2 2(<RP) 3 G) . 

It follows that 

H 2 ' 4 = -2hh^ 1 h 1 h2W 2 wlw 2 2 w 2 3 P i G2 - 2K il) h^ 1 h 1 h 2 w 2 wlw 2 .w 2 P 4 G 1 

O O "i O O 0000/i-> \ (1) 0000 

+2w h-i h l h2Wiw 2 w 3 P4G — h\h2W WiW 2 w 3 (o U2 hi) P4G — 2K_(hh\h2W w 1 w 2 w 3 P^G 
= -2hh-ihih 2 w 2 w\wlwlP A G2 - 2K ( - 1) h- 1 h 1 h 2 w 2 w\wlwlP A Gi 
+2w 2 h-i h^wlwlwlK^ P 4 G - 2K m 1 hh l h 2 w 2 w{wlw 2 3 PiG . 
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Using relation (14), it is easy to see that H 1 ' 4 — H 2 < 4 is symmetric with respect to P and G. In a 
similar way, one can check the symmetric property of the remaining terms. ■ 

4.3 Locality of symmetries generated by recursion operators 

The sufficient conditions for Nijenhuis pseudo-differential operators to produce local objects for- 
mulated in [33, 32, 25] are valid for pseudo-difference operators over the field T s . In a recent paper 
[25], we proved that Nijenhuis operators, which are the product of weakly nonlocal Hamiltonian 
and symplectic operators [18], generate hierarchies of commuting local symmetries and conserved 
densities in involution under some easily verified conditions. To be self-contained, we restate this 
result in [25] for Hamiltonian operator % (26) and symplectic operator (27). 

Consider a Hamiltonian operator T~L of the form (26) and a symplectic operator (27) 
such that = % ■ Z is a Nijenhuis operator (Proposition 3). Assume that 

Then the vector fields p( l >fi = (HZy G rj commute and uA J ) = Zp^'^ G are 
closed 1-forms for i = 1, 2 and j = 0, 1, 2, ■ ■ ■ . 

Theorem 3. For operator 9\ defined by (25), all 9VK^ commute and all 1-forms ZVVK^ are 
closed for i = 1, 2 and j = 0, 1, 2, • • • . 

Proof. We know from (16) that and are commuting generalised symmetries. There- 
fore, we need to verify that L K (i)Z = L K (i)% = while we shall skip the proof of the same 
properties for K^ 2 \ 

We can write anti-symmetric operators Z and H as Z = Y — Y* with operator Y = ^S, and 
H = Z - Z* with 

Z = h- 1 h h lW 2 wl S+2 K (2) A' 1 o K (1) . 

We denote operator D Y [K^] + D* K(1) Y + YD KW by F. It follows from (49) that L kW Z = F-F*. 
Similarly, denoting operator Dz[K^] — D kW Z — ZD* K{1) by G, we have that L K ( 1) 'H = G — G* . 

We now compute operators F and G. The Frechet derivative of and its adjoined operator 
are given by 

D K(1) = -h^w 2 S + d UQ K w + hw 2 S- 1 and D* R(1) = h lW \ S + d UQ K w - h^ 2 w 2 _ x 
Substituting them into the expression for F, we obtain 

F = -h~ 2 (K[ 1] d Ul h + K^d uo h)S + h^d^K^S - br\h- 2 w\ + h^d^K^S + h^h^l 
= —hr\h^2W 2 _i + h _1 hiw 2 , 

where we also used K^> = hw — \d Ul h and = hw\ + \d UQ h from (14). Hence 

L K(1) Z = F - F* = 0. (58) 

By a similar calculation, we find operator G being a symmetric difference operator. It leads to 

L KW H = G-G* = 0. (59) 

Thus we proved the statement. ■ 
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4.4 Yamilov's discretisation of the Krichever-Novikov equation 

It follows from (58) and (59) that 

L KW m = (L KW H)I + HL KW I = . 

This implies that operator 9\ is a recursion operator of differential-difference equation = 
Substituting Q (1) into h defined by (11) and then h into (14) we can rewrite the latter equation 
in the form 

R(m,u,u-i) 

u tl = , (60 

U\ — U-i 

where R is a polynomial 

R(u, v, q) = (av 2 + 2/3v + j)uq + (/3v 2 + Xv + 6) (u + q) + -/v 2 + 2Sv + e 

and 

a = 2{a\ — 0103); /3 = 0205 + 0204 — 0203 — ai<2 6 ; 7 = 2(0405 — 0206); 

A = a\ + a 2 — a| — 0107; 5 = a^aQ + 0506 — 0306 — 0207; e = 2(a@ — 0307), 

where are constant parameters for the Viallet equation (1). Equation (60) can be identified as 
Yamilov's discretization of the Krichever-Novikov equation (YdKN) [9], cf. equation (V4) when 
v = in [10]. Such relations for all the ABS equations and their generalisations introduced in [26] 
was established in [34]. 

It is straightforward to check that R(u\, u, u±) is a symmetric and bi-quadratic polynomial and 
is related to h defined by (11) as follows 

h(u,ui) = ^R(u l} u,ui). 

Thus we can express K( 2 ) (15) , recursion operator 9\ (25), Hamiltonian operator (26) and symplectic 
operator (27) in terms of polynomial R(u,v,q). For example 

Ut2 = = ( _J_ + 1 ) (62) 



(61) 



4 {u\ — U-\) 2 \U2 — U U — U-2 , 

is the next member in the hierarchy of commuting symmetries of the YdKN equation (60), which 
was first given in [28]. Higher symmetries can be obtained by application of the recursion operator 
to the seeds (60) and (62). The locality and commutativity of these symmetries are guaranteed 
by Theorem 3. Cosymmetries of equation (60) coincide with its covariants, i.e. the variational 
derivatives of the conserved densities. Conserved densities can be obtained as residues of the powers 
of the recursion operator. They coincide with the conserved densities for the Viallet equation (1) 
[19]. Obviously, equations (60) and (62) and every member of the hierarchy are multi-Hamiltonian 
systems. For example equation (60) can be written in a Hamiltonian form 

where the symplectic operator is of the form 

2 „ „ 1 2 
-S-S 



R(u\,u,ui) R(ui,u,u\) 

and the Hamiltonian Hq = ln(u\ —u-i) — ^]nR(u,u-i,u) = — ^res ln9^. It follows from Theorem 
3 that Z n = 19V 1 are symplectic operators. The compatibility of these Hamiltonian structures 
follows from the Nijenhuis property of the recursion operator. 
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